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D , Abstract 



Intriguing B — ► ivTjy' decays provide a unique opportunity to study 
a joining of two-gluon configurations arising from the penguin b — * sG 
and b — > sGG transitions, with those inherent to the 77' particle. We 
employ the heavy-light chiral quark model, applied previously to a 
somewhat related B — * Drf decay, as a calculational tool accounting 
for the nonperturbative soft gluon contributions to the amplitude at 
hand. Thereby we arrive at a novel contribution to the singlet pen- 
guin amplitude, which within our model accounts for ~ 10% of the 
measured B — * Kg' amplitude. 
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1 Introduction 



The data on rare decays from CLEO, BaBar, Belle and Tevatron hold promise 
for deepening our understanding of the interplay between flavour- changing 
and QCD dynamics. In particular, the surprise regarding the B — > Krj' am- 
plitude, existing already for considerable time [1, 2], calls for an explanation. 
Namely, the branching ratio for the decay mode B + — > K + r)' is measured 
to be almost six times bigger than the one for B + — > K + tt°, although the 
same basic b — > s penguin mechanism is expected to drive both processes. 
Apparently, this mechanism may have a different appearance when instead of 
the flavour octet pion there is an (almost) flavour singlet rj' particle involved. 

It is well known that some extraordinary properties of the rf particle 
are related to the QCD anomaly. Therefore, a suggestion to explain the 
enhancement of the B — > Krj' amplitude by the QCD anomaly at first sight 
looks very intriguing [3]. 

Our preceding investigation in this direction found that the b — > si]' 
amplitude in the hard gluon regime represents a well defined short distance 
(SD) mechanism [4], but of minor numerical importance for the process at 
hand. The inability of this SD mechanism to account for the measured 
amplitude invites us to explore here the complementary long-distance (LD) 
mechanism. There are conclusions in the literature [5] that the "singlet 
penguin" amplitude (in the language of SU(3) diagrammatic approach [6]) 
contributes substantially to the B — > Krj' enhancement. One of the purposes 
of this note is to investigate this from another angle than done previously, 
i. e. to identify the singlet penguin contribution and estimate it within the 
well-defined microscopic dynamical framework. 

In order to deal with the low energy properties of some of the involved 
gluons and to calculate the physical B — > Krj' amplitude, we have to intro- 
duce an appropriate low-energy description. We will employ ideas based on 
the chiral quark model (xQM)[7, 8] which has been used to describe K — > 2n 
decays [9, 10]. Here we will rely on the extension of such models, namely 
the heavy-light chiral quark model (HLxQM) [11] which has, as the xQM, 
turned out to be a convenient calculational tool for addressing soft-gluon 
contributions, that can be expressed in terms of gluon condensate effects 
[8, 9, 10, 11, 12]. One should note that when the quarks and soft gluons 
in the HLxQM are integrated out, one obtains standard heavy light chiral 
perturbation theory (HLxPT) [13]. The HLxQM has been applied to B — B 
mixing [14] and to decays of the type B — > DD [15]. It has also been applied 
to the decay mode B — > Drf [16] which has similar aspects as the mode we 
consider in this paper. In principle, the HL^QM naturally accounts only 
for soft kaons in the final state. However, as we show in the next section, 



2 



1 K 
i 

i 




Figure 1: A mechanism for B — > Krf due to soft gluons. The crosses corre- 
spond to the gluon condensate. 

assuming the general form of the relevant form factor enables one to perform 
extrapolation from the soft to the hard kaon case. 

2 Penguin contribution to B — >• Krf in the 
heavy-light chiral quark model 

The LD mechanism that we propose for the B — > Krj' decay is shown in 
Figure 1. It accounts for the contribution obtained when a soft gluon (G) 
is emitted from the B — > K transition together with the virtual gluon (G*) 
associated with the penguin b — > s transition. In addition, one of the gluons 
(G 1 ) in the r]'GG vertex in Figure 1 is also assumed to be soft, so that two 
soft gluons form a vacuum condensate. The remaining off-shell gluon from 
the b — > sG* penguin is propagating into the rj'. Now, the B — > KGG* vertex 
denoted by a large circle in Figure 1 can be calculated in the soft K limit 
within the HLxQM, as displayed in Figure 2. 

To begin with it, let us recall that the involved vector current form factors 
in heavy quark physics are defined by 

(K\*fb\B) = U(q 2 )(PB + P KY + f-(q 2 )(PB - pkY 

= F B {q 2 )M B v» + F K {q 2 )p» K) ^ 

where 

F BtK = f + ±f_, and p B = M B v». (2) 

In the M B — * oo limit, dominates, as it can be seen from the scaling 
properties F K ~ \/M B and F B ~ \j\JM B [17]. In addition, considering the 
soft K limit within heavy-light chiral perturbation theory (HLxPT) [13], F K 
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is dominated by the B* pole, and is given by 



(f k ) ft = a ( ^L^ , (3) 



where au = fByMs/ {C 1 + C v ), and C 7 ~ 1, C v ~ are the coefficients 
determined by QCD renormalization of the weak heavy-light current [18]. 
The soft K limit is of course unphysical in our case, and to overcome this we 
employ a double pole structure 

Such a structure, proposed in [19], seems to fit very well the existing data 
and the theoretical requirements on the heavy-light vector form factor, for, 
say, Mi = Mb* , and some parameter 7 fitting M| = 7Mj|* . To determine 7 
we observe that (4) in our limits reads 

Fk{<1 Ut ~ 2(7-1) vp K ' (5) 
so that a comparison with (3) gives 

7^(0) C, f E 



2(7 - 1) C, + C v f w V2 9A ' (6) 

which is general within HLxPT. Knowing the value for F K (0), 7 can be 
determined. We will use the result of the QCD sum rules on the light-cone 
analysis, Fk{0) = 0.34 ± 0.05 [20], implying 7 = 1.27 ± 0.08, in agreement 
with lattice fits [19]. Thus, extrapolating from the soft K to the general case 
we obtain the substitution rule 



1 2(7-1) 1 V2UF K (q 2 ) 



vpk M Bl (i__£.)(i_^y f B9A M B 

B s Bf, 



(7) 



Below we will assume that the form factor for the B — > KGG* vertex also 
has the dipole form (4), because of B* pole dominance in both cases. This 
assumption will hold within HLxPT in the region where it is valid. Therefore 
we will adopt the rule (7) also for the B — > KGG* form factor. For this case, 
the position of the second pole might be somewhat different, which means 
that 7 and F K (0) in (6) should be replaced by 7 G and F^(0), respectively. 

The b — > sG* penguin operator at the quark level is 

g s G P s^Lt a b{DG)l , (8) 
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Figure 2: B — > KGG* in the soft kaon limit. The black box denotes the 
b — > sG* penguin transition. 

where D means a covariant derivative, and 

G r = ^4^KV tb (Fl-F?). (9) 

For the quantity (F{ — F() we take the one loop result [21], 0.26 — ( — 
(2/3) In m 2 c /M^ v \ m —5.2. This result might be slightly changed by pertur- 
bative QCD effects like in [22] for s — > dG*, but we do not enter such details 
here. 

Note that the contribution of the dipole penguin operator F 2 m b sa^ p Rt a bG a ^ u 
is suppressed by the small form factor F 2 ~ 0.2, and we neglect it. 

The bosonization of the coloured quark current in (8) with emission of 
an additional soft gluon is known [14], 

q Ll »Lt a Q v — G H g 8 Glp T±l^>fLH v [A^ + A 2 a a ^ ■ v] j , (10) 

where 

h v = 1 -±^(b; 1 »-iB, 1 ,) (ii) 

is the heavy meson "superfield" [13] and £ = exp(ill/ f n ), with II being the 
standard Goldstone boson 3x3 matrix field. Furthermore, Gh is the meson- 
quark coupling given by G 2 H = 2mp/ f%, where m = 0.250 ± 0.025 GeV is 
the constituent light quark mass, p ~ 1.1 is a hadronic parameter [14] and 
qa = 0.59 ± 0.08 is the axial coupling of Goldstone bosons to heavy mesons. 
Moreover, within the HLxQM we find 
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and I2 is a logarithmically divergent loop integral which is expressed in terms 
of /„. and the gluon condensate [8, 9, 10, 11], as follows 

fl = -^NJ 2 + ^{^). (13) 

It should be noted that the structure in (10) should be rather general, while 
Gh and the explicit expressions for A 12 are model dependent. 

Taking the vector (B*) part of the H v and connecting with the B* prop- 
agator in Figure 2, we obtain in the soft K limit the amplitude 



M(B - KGG*) so { t = G H G P g s G°, (£>(?)• 

x{-(A 1 + A 2 y ua P + 2A 2 v»e Xual3 v x } . (14) 

In order to obtain the general amplitude for B — > KGG* from this equation, 
we perform the substitution (7) for (v ■ px) in the denominator above with 
F K (q 2 ) replaced by F^(q 2 ). 

Concerning the r]'G*G' interaction, it has the general form already used 
in ref. [16], 

V(G* - GV) = -\F, lgg {qV c ^ R ^7' GC ^ > ( 15 ) 

and several groups [23, 24] calculated the form-factor F v i gg (q 2 ) in the pertur- 
bative QCD approach. Since this approach becomes unreliable for gluon mo- 
menta of ~ 1 GeV, we adopt a formula from [24] which interpolates between 
the perturbative QCD region and the anomaly value for zero momentum. 
This formula gives F v > gg (m 2 ,) = 1.55 ± 0.40 GeV -1 , where an error of 25 % 
has been allowed. Accordingly, together with the value of the gluon conden- 
sate, this is the major source of uncertainty in our result below. Taking now 
the vacuum expectation value of the two soft gluons 

g 2 s G a af3 G c XK — > 12(N 2 - 1) \ 9aX9l3K ~ 9aK9f3\) (-^G 2 ) > ( 16 ) 
we obtain the final amplitude 



M(B Kj)m = ^§^^ (^)Ml{A M+A A , 

6JbvMb 7T <- J 

(17) 

where q 2 = m 2 ,. 
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Numerically, with (fG 2 ) = (0.32 ± 0.02 GeV) 4 , f B = 0.18 ± 0.03 GeV, 
and F^(q 2 ) = F K (q 2 ) given by the numbers below Eq. (6), we get 



\M(B -> Krf) m \ = (8 ±3) x 10~ 9 GeV . 



(18) 



It should be noted, that even if the uncertainty of 7 is increased to 50 % when 
replaced by 7g, it has no significant impact on the form factor F^(q 2 ) at the 



that our assumption for the form factor F^(q 2 ) should be rather sound. 

Our result (18) should be compared to the experimental amplitude \M(B - 
-^V)ex P | = (88 ±2) x 10~ 9 GeV. Thus, according to our analysis, only of order 
10% of the B — > Kvj rate enhancement can be ascribed to this gluonic cre- 
ation of 77'. Some additional mechanisms are necessary, such as constructive 
interference of amplitudes for creating rj' in dd and ss state [25, 26]. How- 
ever, the gluonic mechanism studied here, with its distinctive flavour-singlet 
nature, seems to go in the direction of the result of the SU(3) symmetry anal- 
ysis in [5] that shows a substantial singlet penguin contribution to B — > Krj' 
amplitude. 

3 Discussion 

Let us comment here on how our result fits into the existing accounts of 
joining the two-gluon configurations arising from the b — > sG and b — > sGG 
transitions with those inherent to the rj 1 particle, in particular on those relying 
on properties of the rj' particle that are related to the QCD anomaly. The 
attempts to explain some puzzling hadronic weak decays by QCD anomalies 
are well known: the A J = 1/2 enhancement in K — > tttt by the trace anomaly 
[27] and the enhancement of B — > r]'X s decay rates by the axial anomaly 



We entered such study by employing the fact that anomaly permeates all 
distance scales, that enables one to study the role of two-gluon anomalous 
configurations from an extreme SD to a truly LD regime. Our recent study 
[4] shows that there is merely a remnant, the anomaly tail, in the extreme 
SD case. It should be noted that this contribution is obtained as a two quark 
operator for b — > si]' and is very different from the LD contribution presented 
in detail in the previous section. 

In addition, we have also identified some other contributions which we 
have found to be negligible. For example, Fritzsch [28] has suggested that an 
effective interaction of the form 



physical point q 2 = m, 




[3, 28, 29]. 



H eS = aagGpSLbRGnyG 11 ' 



(19) 
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might contribute significantly. We will describe elsewhere [30] different per- 
turbative and nonperturbative contributions to such an effective interaction 
stemming from anomalous two-gluon configurations. Note that already a 
rotation of an appropriate term from Simma and Wyler's paper [31] to 
Fritzsch's form enables one to read off a perturbative contribution to co- 
efficient a above, ~ 4 x 10~ 4 GeV -1 , which gives an amplitude 3-4 
times smaller than (18). 

After observing a systematic suppression of the anomalous two-gluon con- 
tributions through all the distance scales [30], we are focusing in the present 
paper to LD nonperturbative gluon configurations that may be phenomeno- 
logically more relevant. Thereby we are considering the low energy contri- 
bution where a gluon condensate accounts for the emission of soft gluons. A 
priori, our calculation would only be valid in the unphysical case where the 
outgoing kaon is soft. However, one can extrapolate to the physical point 
by introducing a dipole form factor for the B — > KGG* transition (as for 
the standard B — > K transition current). As a result we find a more signifi- 
cant contribution from this mechanism, which can account for ~ 10% of the 
measured amplitude. 

This result has to be compared with findings of [26] that flavour singlet 
contributions to B — > Krj' may be marginal. However, due to quite large un- 
certainties in both amplitudes these two results are actually not inconsistent. 
Note that the major portion of the singlet penguin amplitude in [26] comes 
from the operators corresponding to singlet quark configurations forming rf 
particle. Recently another analysis within SCET appeared [32]. This analysis 
concludes that the "singlet penguin" contribution is essential to understand 
the process B — > Krj' . Unfortunately, a direct comparison between our treat- 
ment of one soft (q 2 ~ 0) and one semi-hard (q 2 ~ m 2 ,) gluon and the one 
by SCET is difficult to perform. Anyway, our contribution corresponding to 
the gluonic configurations forming the gluon condensate is a novel one, and 
may significantly increase the role of the singlet penguin mechanism in the 
direction of the result based on the SU(3) symmetry analysis [5]. 

A number of authors already used the surprising B — > Krj' enhance- 
ment to infer on the contributions beyond the Standard model. However, at 
this stage we need first to consider possible contributions from the specific 
mechanisms within the Standard model, like the one presented here. This 
mechanism seems to provide an additional contribution to the singlet penguin 
topology, the understanding of which may be of importance for explaining 
the data on CP asymmetries in penguin dominated modes [2]. 
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